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Abstract 

We develop a holographic (bottom-up) gravity model for QCD to understand the connection between the peak 
in the trace anomaly and the magnitude of heavy quark energy loss in a strongly-coupled plasma. The potential of 
the scalar field on the gravity side is constructed to reproduce some properties of QCD at finite temperature and its 
parameters are constrained by fitting lattice gauge theory results. The energy loss of heavy quarks (as predicted by 
the holographic model) is found to be strongly sensitive to the medium properties. 
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1. Introduction 

One of the most important properties of the quark-gluon plasma created in heavy ion collisions at the Relativistic 
Heavy Ion Collider is its near-ideal fluid behavior, indicated by the small ratio of shear viscosity to entropy density 
77/s im, which we assume to imply strong coupling. A very useful tool to study strongly coupled gauge theories 
is the AdS/CFT correspondence: in its original form f2\, AdS/CFT correspondence is an equivalence of a (3+1)- 
dimensional N - A super- Yang-Mills (SYM) gauge theory and the type IIB string theory on AdSs x spacetime. 
Using this conjecture and by taking the limit Nc '» A '» I one can study this gauge theory at strong coupling by 
studying classical, two-derivative (super)gravity. 

Because of this very useful aspect of the AdS/CFT correspondence, one is tempted to use thermal N - 4 SYM to 
study and understand properties of finite-temperature QCD. However, there are significant differences between these 
two theories, most of them which originate from the conformal invariance of AT = 4 SYM: this gives rise to a trivial, 
temperature-independent speed of sound (c, - 1/ V3), vanishing bulk viscosity, absence of phase transition, etc. In 
addition, lattice gauge theory results for the QCD trace anomaly |3| show significant violation of the conformal in- 
variance near the crossover temperature Tc- Therefore, in order to capture these important phenomenological features 
of finite-temperature QCD, we are led to consider conformally non-invariant gravity theories with thermodynamics 
properties that can match the relevant QCD phenomenology. 



2. Constructing the potential 

According to the holographic dictionary of the original AdS/CFT correspondence, a scalar field (dilaton) (p in the 
bulk of AdSs space is dual to the TrF^ operator in the TV = 4 SYM at the boundary of that space. Assuming that 
this dictionary holds in general and that the thermodynamics in large-A^c QCD is dominated by the adjoint degrees 
of freedom, we will look for models with a non-trivial dilaton profile, with aim to reproduce the non-trivial QCD 
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thermodynamics. The simplest bottom-up effective reaUzation of a gravity theory dual to a non-conformal gauge 
theory is a 5-dimensional gravity theory coupled to a scalar field (dilaton): 

S = ^^jd'x^^G^R- ^-{dcpf - y(0)J (1) 

where V{(p) is the dilaton potential, is the 5-dimensional gravitational constant and G^y is the 5-dimensional metric. 

The first constraint will come from the requirement that the potential should give asymptotically (0 — » 0) AdS 5 
spacetime, which translates into the conformal invariance of the dual field theory in UV. In order to obtain the AdS 5 
spacetime with radius L as — > 0, the potential must asymptote to a negative constant, so in general it needs to have 
the following form; 

V(4>)^~+^-mY+0(^') (2) 

Such a potential translates into a relevant deformation of the conformal field theory: X.cft — > -Ccft + ^'I '^O^. Here 
is the scale of the deformation and A is the UV dimension of the field theory operator dual to cp. The dimension 
of the operator is determined from A(A -4) = m^L^ with A > 2 (the larger root) and A < 4 in order for the deformation 
to be relevant. This requkement of asymptotically conformal field theory originates from the idea |4| of making the 
connection with QCD by matching the dimension A of to the dimension of TiF^ in QCD at some UV scale Q; 
this UV matching to QCD essentially means that the asymptotic freedom of QCD gets replaced (or approximated) by 
conformal invariance. Therefore, our dual theory will not be asymptotically free, which is also already visible from 
the fact that 77/i is 1 /{An) at all temperatures Q, since we neglected the higher derivative terms in the action ([T]i. This 
in turn means that the validity of our model should not be extended to temperatures too far above or too far below the 
crossover temperature T^. 

A natural choice for one of the terms in the potential is cosh(7'0), since this term interpolates between a constant 
potential for small (p (the CFT limit) and e'"^ (which gives a constant, y-dependent speed of sound) for large (p and, as 
discussed in |6|, should give a temperature -dependent speed of sound. The coefficient y gives the value of the speed 
of sound in the IR (and will be numerically determined by fitting the lattice results), but it also fixes the dimension A, 
according to (j2]l. In order to have more freedom in choosing this dimension, we will introduce an additional <p^ term 
in the potential. We will also introduce (p* and (p^ terms to fine-tune the fit of the cross-over behavior of the lattice 
results. Finally, our potential has the following form ||4]|6l: 

V(cf>) = ;^ (-12 cosh(r0) + ^20' + b4<p' + be4>^) . (3) 

A similar bottom-up approach has been studied in |7ji8|, where the authors explored the role played by the gauge 
theory's beta function in the construction of the scalar field potential in the gravity theory. 



3. Fitting tlie lattice 

We will use the following ansatz for the metric in the Einstein frame 

ds^ = G.yd^dx' = e2A(r) (_/,(^)^,2 ^ ^^2^ ^ ^2B(r) ^ ^ ^ ^^^^^ 

^ ■' h{r) 

which is dictated by the required symmetries |6i and where the SO(3,l) boost invariance in (t,x) - {t,x\,X2,Xi) 
direction is broken by the finite temperature, i.e. presence of the black hole with regular horizon at some finite r - rn 
defined by /i(r/,) - 0. In this ansatz we still have some gauge freedom left to reparametrize the radial direction r, so 
we will follow ^ and use the gauge choice 0(r) = r, which means that the conformal (UV) boundary is at r = 0. 
With this gauge choice and ansatz Q we can solve the equations of motion from the action ([T]i with potential of the 
form ([3]l using the method developed in |6| in which one constructs a nonlinear 'master' differential equation for a 
generating function G{4>), from which one then obtains all three unknown functions in Q, A{(p), B{(p) and h{(p). 

Once we solve the equations of motion, we can find the entropy density and temperature using Hawking's formu- 
las: 

,^^,iMn)^ T ^^e^('"^-''^'-"Ah'{rH)\. (5) 

Kt 4-7T ' ' 
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From here we can then find the speed of sound as Cj - dlog T/dlog s and other thermodynamic quantities. Since 
we now have the speed of sound cj(r) dependent on the parameters of the potential we can numerically deter- 
mine these parameters by fitting our results to the lattice QCD data (with dynamical fermions) from the hotQCD 
collaboration jSj. Satisfactory fits were obtained by using the following values for the parameters in the potential: 
y = 0.606, b2 - 0.703, ^4 = -0.12, be = 0.00325, which corresponds to dimension A a; 3. Our fit together with 
the lattice results is shown in Figure 1 (leftmost plot) and after fitting the speed of sound all other thermodynamic 
quantities are automatically fitted as well (e.g. trace anomaly, center plot). In the plots also indicated is the critical 
temperature Tc, which we define as the temperature where the speed of sound has a minimum. One must also define a 
rule for converting the temperature in GeV (as a function of which the lattice results are reported) and the temperature 
in 1/L units, or the 'AdS' units (as a function of which one obtains holographic results), in order to be able to compare 
both of them on the same plot. We do that by finding r„„„ in AdS units, where the holographic speed of sound reaches 
minimum, and by postulating a linear conversion rule T{GeV)/Tc = T(AdS)/T,„i„, we obtain the same shapes of the 
thermodynamical quantities as functions of temperature as the ones from the lattice. 




T I GeV T I GeV T / GeV 

Figure 1. Comparison of the speed of sound (left), trace anomaly 8 = e - 3p (center) and the renormalized Polyakov loop (right) as a function 
of temperature from our non-conformal holographic model (red curves) and the lattice results from the hotQCD collaboration |3| (dashed curves 
and points). Also indicated is the critical temperature which we define as the temperature where the speed of sound reaches minimum. 



4. Heavy quark energy loss 

We now introduce a quark of mass niQ, which is dual to a curved string in the bulk stretching from a D4 brane at 
~ l/niQ to the horizon rfj |9|. We are working in the probe approximation, where one neglects the backreaction 

of the metric due to the introduction of the string. Assuming quantum corrections to be negligible, the dynamics of 

the string is governed by the classical Nambu-Goto action: 

where a' = f^, squared fundamental string length, cr" are the string worldsheet coordinates, X'' are its spacetime 
coordinates and q{(p) is the coupling between the string and the dilaton 0, which originates when one expresses the 
string frame metric (which directly enters the Nambu-Goto action) in terms of the Einstein frame metric (j4]l. For this 

coupling we use the result from the five-dimensional non-critical string theory q{cf)) - e [7J. According to the 
gauge/gravity dictionary, Polyakov loops are related to the string action: 

\{L(T))\ = e-^e(TVT ^ ^-S,a(0) (7) 

and since the Nambu-Goto action ^ contains a' (at this moment the only undetermined parameter in our model), we 
can use the lattice results for the Polyakov loop [3 ] and by fitting the prediction of our model to the data (using the 
procedure from |.10i ), determine a' . In that way, our model is completely constrained by the lattice data and the energy 
loss then becomes a prediction. From the fit in the rightmost plot in Figure 1, we see that we get roughly L? ja' ~ 1. 
This means that, in principle, c' -corrections (in the string worldsheet fluctuations and in the eff'ective gravity action) 
can be important in this case, and they must be explicitly calculated in order to see their elTect on thermodynamics 
and on the energy loss. 
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Assuming that the quark is heavy enough and that the range of relevant temperatures is low enough so that r„, <K 
rn, we can evaluate the quark energy loss by using the traihng string model 1(121 : 

(8) 

ax Zna 

where v is the velocity of the quark (assumed uniformly moving) and r» is defined as /i(r») = v^. Of course, in the 
CFT limit this expression becomes the famihar drag force ~ ^/AT^yv ifTZl . 

Assuming that the trailing string model is applicable to both charm and bottom quarks (in the sense that r,„ <sc rn) 
and that they obey the usual relativistic dispersion relation p - ymv, we can compute their energy loss in our non- 
conformal holographic model using the formula ([8]l. Figure 2 shows the energy loss as a function of temperature for a 
fixed momentum and as a function of momentum for several fixed temperatures. In all plots the energy loss from the 
non-conformal model is compared to the CFT limit of our model, where the same ja' ~ 1 was used. 




T I GeV T I GeV p I GeV 

Figure 2. Energy loss for charm and bottom quarks (and their ratios) as a function of the temperature for fixed momentum p = 10 GeV (left and 
center) and as a function of momentum for several fixed temperatures (right) predicted in our non-conformal holographic model, compared to 
the CFT limit of the model. 



From Figure 2 we see that a significant deviation from the CFT results is found for temperatures around Tc, where 
the trace anomaly is maximal. This difference should affect the relative ratio of the energy loss observables, such 
as R/^A, since we expect that a quark traveling through a realistic, expanding quark-gluon plasma spends significant 
amount of time in the temperature region around Tc- 

5. Conclusions and outlook 

Gauge/string duality gives us means to address the violation of conformal invariance in a strongly coupled plasma 
near Tc and see how the trace anomaly and the details of equation of state affect the energy loss of heavy quarks. We 
showed that, in principle, the details of equation of state are important for the heavy quark energy loss in the crossover 
region (near Tc) and may affect observables such as the nuclear modification factor of jets in heavy ion collisions (\3\. 

We thank W. Zajc and A. Dumitru for helpful discussions. We acknowledge support by US-DOE Nuclear Science 
Grant No. DE-FG02-93ER40764. 
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